Fine mapping of TFL, a major gene regulating fruit length in snake gourd (Trichosanthes anguina L).

Fruit length is a crucial agronomic trait of snake gourd (Trichosanthes anguina L); however, genes associated with fruit length have not been characterised. In this study, F2 snake gourd populations were generated by crossing the inbred lines, S1 and S2 (fruit lengths: 110 and 20 cm, respectively). Subsequently, bulk segregant analysis, sequencing, and fine-mapping were performed on the F2 population to identify target genes. Our findings suggest that the fruit length of snake gourd is regulated by a major-effect regulatory gene. Mining of genes regulating fruit length in snake gourd to provide a basis for subsequent selection and breeding of new varieties. Genotype-phenotype association analysis was performed on the segregating F2 population comprising 6,000 plants; the results indicate that the target gene is located on Chr4 (61,846,126-61,865,087 bp, 18.9-kb interval), which only carries the annotated candidate gene, Tan0010544 (designated TFL). TFL belongs to the MADS-box family, one of the largest transcription factor families. Sequence analysis revealed a non-synonymous mutation of base C to G at position 202 in the coding sequence of TFL, resulting in the substitution of amino acid Gln to Glu at position 68 in the protein sequence. Subsequently, an InDel marker was developed to aid the marker-assisted selection of TFL. The TFL in the expression parents within the same period was analysed using quantitative real-time PCR; the TFL expression was significantly higher in short fruits than long fruits. Therefore, TFL can be a candidate gene for determining the fruit length in snake gourd. Collectively, these findings improve our understanding of the genetic components associated with fruit length in snake gourds, which could aid the development of enhanced breeding strategies for plant species.

The microRNA miR482 regulates NBS-LRR genes in response to ALT1 infection in apple.

Plants are frequently attacked by a variety of pathogens and thus have evolved a series of defense mechanisms, one important mechanism is resistance gene (R gene)-mediated disease resistance, but its expression is tightly regulated. NBS-LRR genes are the largest gene family of R genes. microRNAs (miRNAs) target to a number of NBS-LRR genes and trigger the production of phased small interfering RNAs (phasiRNAs) from these transcripts. phasiRNAs cis or trans regulate NBS-LRR genes, which can result in the repression of R gene expression. In this study, we screened for upregulated miR482 in the susceptible apple cultivar 'Golden Delicious' (GD) after inoculation with the fungal pathogen Alternaria alternata f. sp. mali (ALT1). Additionally, through combined degradome sequencing, we identified a gene targeted by miR482, named MdTNL1, a gene encoding a TIR-NBS-LRR (Toll/interleukin1 receptor-nucleotide binding site-leucine-rich repeat) protein. This gene exhibited a significant down-regulation post ALT1 inoculation, suggesting an impact on gene expression mediated by miRNA regulation. miR482 could cleave MdTNL1 and generate phasiRNAs at the cleavage site. We found that overexpression of miR482 inhibited the expression of MdTNL1 and thus reduced the disease resistance of GD, while silencing of miR482 increased the expression of MdTNL1 and thus improved the disease resistance of GD. This work elucidates key mechanisms underlying the immune response to Alternaria infection in apple. Identification of the resistance genes involved will enable molecular breeding for prevention and control of Alternaria leaf spot disease in this important fruit crop.

Identification and expression analyses of B3 genes reveal lineage-specific evolution and potential roles of REM genes in pepper.

BACKGROUND: The B3 gene family, one of the largest plant-specific transcription factors, plays important roles in plant growth, seed development, and hormones. However, the B3 gene family, especially the REM subfamily, has not been systematically and functionally studied. RESULTS: In this study, we performed genome-wide re-annotation of B3 genes in five Solanaceae plants, Arabidopsis thaliana, and Oryza sativa, and finally predicted 1,039 B3 genes, including 231 (22.2%) newly annotated genes. We found a striking abundance of REM genes in pepper species (Capsicum annuum, Capsicum baccatum, and Capsicum chinense). Comparative motif analysis revealed that REM and other subfamilies (ABI3/VP1, ARF, RAV, and HSI) consist of different amino acids. We verified that the large number of REM genes in pepper were included in the specific subgroup (G8) through the phylogenetic analysis. Chromosome location and evolutionary analyses suggested that the G8 subgroup genes evolved mainly via a pepper-specific recent tandem duplication on chromosomes 1 and 3 after speciation between pepper and other Solanaceae. RNA-seq analyses suggested the potential functions of REM genes under salt, heat, cold, and mannitol stress conditions in pepper (C. annuum). CONCLUSIONS: Our study provides evolutionary and functional insights into the REM gene family in pepper.

Contiguous identity between entire coding regions of transgenic and native genes rather than special regions is essential for a strong co-suppression.

The discovery of co-suppression in plants has greatly boosted the study of gene silencing mechanisms, but its triggering mechanism has remained a mystery. In this study, we explored its possible trigger mechanism by using Fatty acid desaturase 2 (FAD2) and Fatty acid elongase 1 (FAE1) strong co-suppression systems. Analysis of small RNAs in FAD2 co-suppression lines showed that siRNAs distributed throughout the coding region of FAD2 with an accumulated peak. However, mutations of the peak siRNA-matched site and siRNA derived site had not alleviated the co-suppression of its transgenic lines. Synthetic FAD2 (AtFAD2sm), which has synonymous mutations in the entire coding region, failed to trigger any co-suppression. Furthermore, 5' and 3' portions of AtFAD2 and AtFAD2sm were swapped to form two hybrid genes, AtFAD2-3sm and AtFAD2-5sm. 80 % and 92 % of their transgenic lines exhibited co-suppression, respectively. Finally, FAE1s with different degrees of the continuous sequence identity compared with AtFAE1 were tested in their Arabidopsis transgenic lines, and the results showed the co-suppression frequency was reduced as their continuous sequence identity stepped down. This work suggests that contiguous identity between the entire coding regions of transgenic and native genes rather than a special region is essential for a strong co-suppression.

Profiling of the gene expression and alternative splicing landscapes of Eucalyptus grandis.

Eucalyptus is a widely planted hardwood tree species due to its fast growth, superior wood properties and adaptability. However, the post-transcriptional regulatory mechanisms controlling tissue development and stress responses in Eucalyptus remain poorly understood. In this study, we performed a comprehensive analysis of the gene expression profile and the alternative splicing (AS) landscape of E. grandis using strand-specific RNA-Seq, which encompassed 201 libraries including different organs, developmental stages, and environmental stresses. We identified 10 416 genes (33.49%) that underwent AS, and numerous differentially expressed and/or differential AS genes involved in critical biological processes, such as primary-to-secondary growth transition of stems, adventitious root formation, aging and responses to phosphorus- or boron-deficiency. Co-expression analysis of AS events and gene expression patterns highlighted the potential upstream regulatory role of AS events in multiple processes. Additionally, we highlighted the lignin biosynthetic pathway to showcase the potential regulatory functions of AS events in the KNAT3 and IRL3 genes within this pathway. Our high-quality expression atlas and AS landscape serve as valuable resources for unravelling the genetic control of woody plant development, long-term adaptation, and understanding transcriptional diversity in Eucalyptus. Researchers can conveniently access these resources through the interactive ePlant browser (https://bar.utoronto.ca/eplant_eucalyptus).

GIGANTEA accelerates wheat heading time through gene interactions converging on FLOWERING LOCUS T1.

Precise regulation of flowering time is critical for cereal crops to synchronize reproductive development with optimum environmental conditions, thereby maximizing grain yield. The plant-specific gene GIGANTEA (GI) plays an important role in the control of flowering time, with additional functions on the circadian clock and plant stress responses. In this study, we show that GI loss-of-function mutants in a photoperiod-sensitive tetraploid wheat background exhibit significant delays in heading time under both long-day (LD) and short-day photoperiods, with stronger effects under LD. However, this interaction between GI and photoperiod is no longer observed in isogenic lines carrying either a photoperiod-insensitive allele in the PHOTOPERIOD1 (PPD1) gene or a loss-of-function allele in EARLY FLOWERING 3 (ELF3), a known repressor of PPD1. These results suggest that the normal circadian regulation of PPD1 is required for the differential effect of GI on heading time in different photoperiods. Using crosses between mutant or transgenic plants of GI and those of critical genes in the flowering regulation pathway, we show that GI accelerates wheat heading time by promoting FLOWERING LOCUS T1 (FT1) expression via interactions with ELF3, VERNALIZATION 2 (VRN2), CONSTANS (CO), and the age-dependent microRNA172-APETALA2 (AP2) pathway, at both transcriptional and protein levels. Our study reveals conserved GI mechanisms between wheat and Arabidopsis but also identifies specific interactions of GI with the distinctive photoperiod and vernalization pathways of the temperate grasses. These results provide valuable knowledge for modulating wheat heading time and engineering new varieties better adapted to a changing environment.

Identification and function analysis of GABA branch three gene families in the cotton related to abiotic stresses.

γ -aminobutyric acid (GABA) is closely related to the growth, development and stress resistance of plants. Combined with the previous study of GABA to promote the cotton against abiotic stresses, the characteristics and expression patterns of GABA branch gene family laid the foundation for further explaining its role in cotton stress mechanism. Members of GAD, GAB-T and SSADH (three gene families of GABA branch) were identified from the Gossypium hirsutum, Gossypium barbadense, Gossypium arboreum and Gossypium raimondii genome. The GABA branch genes were 10 GAD genes, 4 GABA-T genes and 2 SSADH genes. The promoter sequences of genes mainly contains response-related elements such as light, hormone and environment.Phylogenetic analysis shows that GAD indicating that even in the same species, the homologous sequences in the family. The GABA-T gene of each cotton genus was in sum the family had gene loss in the process of dicotyledon evolution. SSADH families Gossypium hirsutum, Gossypium barbadense, Gossypium arboreum and Gossypium raimondii were closely related to the dicot plants.GABA gene is involved in the regulation of salt stress and high temperature in Gossypium hirsutum.GABA attenuated part of the abiotic stress damage by increasing leaf protective enzyme activity and reducing reactive oxygen species production.This lays the foundation for a thorough analysis of the mechanism of GABA in cotton stress resistance.

Dissection of a rapidly evolving wheat resistance gene cluster by long-read genome sequencing accelerated the cloning of Pm69.

Gene cloning in repeat-rich polyploid genomes remains challenging. Here, we describe a strategy for overcoming major bottlenecks in cloning of the powdery mildew resistance gene (R-gene) Pm69 derived from tetraploid wild emmer wheat. A conventional positional cloning approach was not effective owing to suppressed recombination. Chromosome sorting was compromised by insufficient purity. A Pm69 physical map, constructed by assembling Oxford Nanopore Technology (ONT) long-read genome sequences, revealed a rapidly evolving nucleotide-binding leucine-rich repeat (NLR) R-gene cluster with structural variations. A single candidate NLR was identified by anchoring RNA sequencing reads from susceptible mutants to ONT contigs and was validated by virus-induced gene silencing. Pm69 is likely a newly evolved NLR and was discovered in only one location across the wild emmer wheat distribution range in Israel. Pm69 was successfully introgressed into cultivated wheat, and a diagnostic molecular marker was used to accelerate its deployment and pyramiding with other R-genes.

GWAS identifies two important genes involved in Chinese chestnut weight and leaf length regulation.

There are many factors that affect the yield of Chinese chestnut (Castanea mollissima), with single nut weight (SNW) being one of the most important. Leaf length is also related to Chinese chestnut yield. However, the genetic architecture and gene function associated with Chinese chestnut nut yield have not been fully explored. In this study, we performed genotyping by sequencing 151 Chinese chestnut cultivars, followed by a genome-wide association study (GWAS) on six horticultural traits. First, we analyzed the phylogeny of the Chinese chestnut and found that the Chinese chestnut cultivars divided into two ecotypes, a northern and southern cultivar group. Differences between the cultivated populations were found in the pathways of plant growth and adaptation to the environment. In the selected regions, we also found interesting tandemly arrayed genes that may influence Chinese chestnut traits and environmental adaptability. To further investigate which horticultural traits were selected, we performed a GWAS using six horticultural traits from 151 cultivars. Forty-five loci that strongly associated with horticultural traits were identified, and six genes highly associated with these traits were screened. In addition, a candidate gene associated with SNW, APETALA2 (CmAP2), and another candidate gene associated with leaf length (LL), CRYPTOCHROME INTERACTING BASIC HELIX-LOOP-HELIX 1 (CmCIB1), were verified in Chinese chestnut and Arabidopsis (Arabidopsis thaliana). Our results showed that CmAP2 affected SNW by negatively regulating cell size. CmCIB1 regulated the elongation of new shoots and leaves by inducing cell elongation, potentially affecting photosynthesis. This study provided valuable information and insights for Chinese chestnut breeding research.

Overexpression of SlCRF6 in tomato inhibits leaf development and affects plant morphology.

Cytokinin response factors (CRFs) are transcription factors (TFs) that are specific to plants and have diverse functions in plant growth and stress responses. However, the precise roles of CRFs in regulating tomato plant architecture and leaf development have not been comprehensively investigated. Here, we identified a novel CRF, SlCRF6, which is involved in the regulation of plant growth via the gibberellin (GA) signaling pathway. SlCRF6-overexpressing (SlCRF6-OE) plants displayed pleiotropic phenotypic changes, including reduced internode length and leaf size, which caused dwarfism in tomato plants. This dwarfism could be alleviated by application of exogenous GA3. Remarkably, quantitative real-time PCR (qRTPCR), a dual luciferase reporter assay and a yeast one-hybrid (Y1H) assay revealed that SlCRF6 promoted the expression of SlDELLA (a GA signal transduction inhibitor) in vivo. Furthermore, transgenic plants displayed variegated leaves and diminished chlorophyll content, resulting in decreased photosynthetic efficiency and less starch than in wild-type (WT) plants. The results of transient expression assays and Y1H assays indicated that SlCRF6 suppressed the expression of SlPHAN (leaf morphology-related gene). Collectively, these findings suggest that SlCRF6 plays a crucial role in regulating tomato plant morphology, leaf development, and the accumulation of photosynthetic products.

Disruption of the rice ALS1 localized in chloroplast causes seedling-lethal albino phenotype.

Chloroplasts are the organelles responsible for photosynthesis and regulate normal plant growth. Although translation elongation factors play important roles in chloroplast development, functional studies of chloroplast translation elongation factors in higher plants remain very sparse. Here, we obtained a rice mutant exhibiting seedling-lethal albino phenotype and named it albino and lethal seedling 1 (als1). Consistently, low content of photosynthetic pigments, malformed chloroplasts and defective photosynthesis were observed in als1 mutant leaves. Map-based cloning experiment showed that als1 mutant had a T base insertion in Os02g0595700, causing a frame shift and premature stop codon. ALS1 encoded a GTP-binding protein EF-Tu, which acts as a translation elongation factor in chloroplast protein translation. ALS1 was found to be expressed throughout plant with highest expression level in young leaves. Moreover, ALS1 was located in chloroplast, whereas the truncated als1 could not normally be located in chloroplast. Additionally, the ALS1 mutation significantly influenced the expression of downstream genes, such as genes relevant to chlorophyll biosynthesis, photosynthesis as well as chloroplast development. These results show that ALS1 acts as a key regulator of chloroplast development and plant growth.

RNA-Sequencing Analysis Revealed Genes Associated with Sweet Potato (Ipomoea batatas (L.) Lam.) Responses to Stem Rot during Different Infection Stages.

The sweet potato, which is an important tuber crop in China, is susceptible to a variety of pathogens and insect pests during cultivation and production. Stem rot is a common sweet potato disease that seriously affects tuber yield and quality. Unfortunately, there have been relatively few studies on the mechanism mediating the stem rot resistance of sweet potatoes. In this study, a transcriptome sequencing analysis was completed using Xushu 48 samples at different stages (T1, T2, and T3) of the stem rot infection. The T1 vs. T2, T1 vs. T3, and T2 vs. T3 comparisons detected 44,839, 81,436, and 61,932 differentially expressed genes (DEGs), respectively. The DEGs encoded proteins primarily involved in alanine, aspartate, and glutamate metabolism (ko00250), carbon fixation in photosynthetic organisms (ko00710), and amino sugar and nucleotide sugar metabolism (ko00520). Furthermore, some candidate genes induced by phytopathogen infections were identified, including gene-encoding receptor-like protein kinases (RLK5 and RLK7), an LRR receptor-like serine/threonine protein kinase (SERK1), and transcription factors (bHLH137, ERF9, MYB73, and NAC053). The results of this study provide genetic insights that are relevant to future explorations of sweet potato stem rot resistance, while also providing the theoretical basis for breeding sweet potato varieties that are resistant to stem rot and other diseases.

Integrated analysis of HSP20 genes in the developing flesh of peach: identification, expression profiling, and subcellular localization.

BACKGROUND: Plant HSP20s are not only synthesized in response to heat stress but are also involved in plant biotic and abiotic stress resistance, normal metabolism, development, differentiation, survival, ripening, and death. Thus, HSP20 family genes play very important and diverse roles in plants. To our knowledge, HSP20 family genes in peach have not yet been characterized in detail, and little is known about their possible function in the development of red flesh in peach. RESULTS: In total, 44 PpHSP20 members were identified in the peach genome in this study. Forty-four PpHSP20s were classified into 10 subfamilies, CI, CII, CIII, CV, CVI, CVII, MII, CP, ER, and Po, containing 18, 2, 2, 10, 5, 1, 1, 2, 1, and 2 proteins, respectively. Among the 44 PpHSP20 genes, 6, 4, 4, 3, 7, 11, 5, and 4 PpHSP20 genes were located on chromosomes 1 to 8, respectively. In particular, approximately 15 PpHSP20 genes were located at both termini or one terminus of each chromosome. A total of 15 tandem PpHSP20 genes were found in the peach genome, which belonged to five tandemly duplicated groups. Overall, among the three cultivars, the number of PpHSP20 genes with higher expression levels in red flesh was greater than that in yellow or white flesh. The expression profiling for most of the PpHSP20 genes in the red-fleshed 'BJ' was higher overall at the S3 stage than at the S2, S4-1, and S4-2 stages, with the S3 stage being a very important period of transformation from a white color to the gradual anthocyanin accumulation in the flesh of this cultivar. The subcellular localizations of 16 out of 19 selected PpHSP20 proteins were in accordance with the corresponding subfamily classification and naming. Additionally, to our knowledge, Prupe.3G034800.1 is the first HSP20 found in plants that has the dual targets of both the endoplasmic reticulum and nucleus. CONCLUSIONS: This study provides a comprehensive understanding of PpHSP20s, lays a foundation for future analyses of the unknown function of PpHSP20 family genes in red-fleshed peach fruit and advances our understanding of plant HSP20 genes.

Why old duplicated genes are not thrown away in (paleo)polyploids? Example from the petC gene in Brassica napus.

Duplication of genes at different time period, through recurrent and frequent polyploidization events, have played a major role in plant evolution, adaptation and diversification. Interestingly, some of the ancestral duplicated genes (referred as paleologs), have been maintained for millions of years, and there is still a poor knowledge of the reasons of their retention, especially when testing the phenotypic effect of individual copies by using functional genetic approaches. To fill this gap, we performed functional genetic (CRISPR-Cas9), physiological, transcriptomic and evolutionary studies to finely investigate this open question, taking the example of the petC gene (involved in cytochrome b6/f and thus impacting photosynthesis) that is present in four paleologous copies in the oilseed crop Brassica napus. RNA-Seq and selective pressure analyses suggested that all paleologous copies conserved the same function and that they were all highly transcribed. Thereafter, the Knock Out (K.O.) of one, several or all petC copies highlighted that all paleologous copies have to be K.O. to suppress the gene function. In addition, we could determine that phenotypic effects in single and double mutants could only be deciphered in high light conditions. Interestingly, we did not detect any significant differences between single mutants K.O. for either the A03 or A09 copy (despite being differentially transcribed), or even between mutants for a single or two petC copies. Altogether, this work revealed that petC paleologs have retained their ancestral function and that the retention of these copies is explained by their compensatory role, especially in optimal environmental conditions.

Hidden prevalence of deletion-inversion bi-alleles in CRISPR-mediated deletions of tandemly arrayed genes in plants.

Tandemly arrayed genes (TAGs) with functional redundancy and chromosomal linkage constitute 14 ~ 35% in sequenced plant genomes. The multiplex CRISPR system is the tool of choice for creating targeted TAG deletions. Here, we show that up to ~80% of CRISPR-mediated TAG knockout alleles in Arabidopsis and rice are deletion-inversion (delinver) bi-alleles, which are easily misidentified as homozygous deletion alleles by routine PCR-based genotyping. This can lead to misinterpretation of experimental data and production of progenies with genetic heterogeneity in an unnoticed manner. In ~2,650 transgenic events, delinver mutation frequencies are predominantly correlated with deletion frequencies but unrelated to chromosomal locations or deletion sizes. Delinver mutations also occur frequently at genomic non-TAG loci during multiplexed CRISPR editing. Our work raises the alarm about delinver mutations as common unwanted products of targeted TAG deletions in plants and helps prevent false interpretation of plant TAG functions due to this hidden genotype issue.

Molecular characterization of stable QTL and putative candidate genes for grain zinc and iron concentrations in two related wheat populations.

KEY MESSAGE: Major QTL for grain zinc and iron concentrations were identified on the long arm of chromosomes 2D and 6D. Gene-based KASP markers were developed for putative candidate genes TaIPK1-2D and TaNAS10-6D. Micronutrient malnutrition is one of the most common public health problems in the world. Biofortification, the most attractive and sustainable solution to surmount malnutrition requires the development of micronutrient enriched new crop cultivars. In this study, two recombinant inbred line (RIL) populations, ZM175/XY60 and ZM175/LX987, were used to identify QTL for grain zinc concentration (GZnC), grain iron concentration (GFeC) and thousand grain weight (TGW). Eight QTL for GZnC, six QTL for GFeC and five QTL for TGW were detected. Three QTL on chromosomes 2DL and 4BS and chromosome 6A showed pleiotropic effects on all three traits. The 4BS and 6A QTL also increased plant height and might be Rht-B1a and Rht25a, respectively. The 2DL locus within a suppressed recombination region was identified in both RIL populations and the favorable allele simultaneously increasing GZnC, GFeC and TGW was contributed by XY60 and LX987. A QTL on chromosome 6DL associated only with GZnC was detected in ZM175/XY60 and was validated in JD8/AK58 RILs using kompetitive allele-specific PCR (KASP) marker K_AX-110119937. Both the 2DL and 6DL QTL were new loci for GZnC. Based on gene annotations, sequence variations and expression profiles, the phytic acid biosynthesis gene TaIPK1-2D and nicotianamine synthase gene TaNAS10-6D were predicted as candidate genes. Their gene-based KASP markers were developed and validated in a cultivar panel of 343 wheat accessions. This study investigated the genetic basis of GZnC and GFeC and provided valuable candidate genes and markers for breeding Zn- and Fe-enriched wheat.

Identification and adaptive evolution analysis of glutaredoxin genes in Populus spp.

Glutaredoxin (GRX) is a class of small redox proteins widely involved in cellular redox homeostasis and the regulation of various cellular processes. The role of GRX gene in the differentiation of Populus spp. is rarely reported. We compared the similarities and differences of GRX genes among four sections of poplar using bioinformatics, corrected the annotations of some GRX genes, and focused on analysing their transcript profiling and adaptive evolution in Populus spp. A total of 219 GRX genes were identified in four sections of poplar, among which annotations for 13 genes were corrected. Differences in GRX genes were found between sect. Turanga, represented by P. euphratica, and other poplar sections. Most notably, P. euphratica had the smallest number of duplication events for GRX genes (n = 9) and no tandem duplications, whereas there were >25 duplication events for all other poplars. Furthermore, we detected 18 pairs of GRX genes under positive selection pressure in various sections of poplar, and identified two groups of GRX genes in the Salicaceae that potentially underwent positive selection. Expression profiling results showed that the PtrGRX34 and its orthologous genes were upregulated under stress treatments. In summary, the GRX gene family underwent expansion during poplar differentiation, and some genes underwent rapid evolution during this process, which may be beneficial for Populus spp. to adapt to environmental changes. This study may provide more insights into the molecular mechanisms of Populus spp. adaptation to environmental changes and the adaptive evolution of GRX genes.

Origin, evolution, and diversification of the wall-associated kinase gene family in plants.

KEY MESSAGE: The study of the origin, evolution, and diversification of the wall-associated kinase gene family in plants facilitates their functional investigations in the future. Wall-associated kinases (WAKs) make up one subfamily of receptor-like kinases (RLKs), and function directly in plant cell elongation and responses to biotic and abiotic stresses. The biological functions of WAKs have been extensively characterized in angiosperms; however, the origin and evolutionary history of the WAK family in green plants remain unclear. Here, we performed a comprehensive analysis of the WAK family to reveal its origin, evolution, and diversification in green plants. In total, 1061 WAK genes were identified in 37 species from unicellular algae to multicellular plants, and the results showed that WAK genes probably originated before bryophyte differentiation and were widely distributed in land plants, especially angiosperms. The phylogeny indicated that the land plant WAKs gave rise to five clades and underwent lineage-specific expansion after species differentiation. Cis-acting elements and expression patterns analyses of WAK genes in Arabidopsis and rice demonstrated the functional diversity of WAK genes in these two species. Many gene gains and losses have occurred in angiosperms, leading to an increase in the number of gene copies. The evolutionary trajectory of the WAK family during polyploidization was uncovered using Gossypium species. Our results provide insights into the evolution of WAK genes in green plants, facilitating their functional investigations in the future.

Refinement of Rf1-gene localization and development of the new molecular markers for fertility restoration in sunflower.

BACKGROUND: Ability to restore male fertility is important trait for sunflower breeding. The most commonly used fertility restoration gene in the production of sunflower hybrids is Rf1. The localization of Rf1 on the linkage group 13 has been previously shown, however, its exact position, its sequence and molecular mechanism for fertility restoration remain unknown. Therefore, several markers linked to Rf1 gene, commonly used for MAS, don't always allow to identify the genotype of plants. For this reason, the search for new markers and precise localization of the Rf1 gene is an urgent task. METHODS AND RESULTS: Based on previously identified single nucleotide polymorphisms (SNPs) at LG13, significantly associated with the ability to restore male fertility, two markers have been developed that have performed well after careful evaluation. These markers, together with other Rf1 markers, were applied for genotyping 72 diversity panel accessions and 291 individuals of F2 segregating population, obtained from crossing the cytoplasmic male sterility (CMS) AHO33 and restorer RT085HO lines. The analysis revealed no recombinants between Rf1 gene and SRF833 marker, the distance between Rf1 and SRF122 marker was 1.0 cM. CONCLUSIONS: Data obtained made it possible to specify the localization of the Rf1 gene and reduce the list of candidate genes to the 3 closely linked PPR-genes spanning a total of 59 Kb. However, it cannot be ruled out that analysis of the candidate region in the genome of fertility restorer lines can reveal new candidate genes in this locus that are absent in the cytoplasmic male sterility maintainer reference sequence.

Comparative pangenome analyses provide insights into the evolution of Brassica rapa resistance gene analogues (RGAs).

Brassica rapa is grown worldwide as economically important vegetable and oilseed crop. However, its production is challenged by yield-limiting pathogens. The sustainable control of these pathogens mainly relies on the deployment of genetic resistance primarily driven by resistance gene analogues (RGAs). While several studies have identified RGAs in B. rapa, these were mainly based on a single genome reference and do not represent the full range of RGA diversity in B. rapa. In this study, we utilized the B. rapa pangenome, constructed from 71 lines encompassing 12 morphotypes, to describe a comprehensive repertoire of RGAs in B. rapa. We show that 309 RGAs were affected by presence-absence variation (PAV) and 223 RGAs were missing from the reference genome. The transmembrane leucine-rich repeat (TM-LRR) RGA class had more core gene types than variable genes, while the opposite was observed for nucleotide-binding site leucine-rich repeats (NLRs). Comparative analysis with the B. napus pangenome revealed significant RGA conservation (93%) between the two species. We identified 138 candidate RGAs located within known B. rapa disease resistance QTL, of which the majority were under negative selection. Using blackleg gene homologues, we demonstrated how these genes in B. napus were derived from B. rapa. This further clarifies the genetic relationship of these loci, which may be useful in narrowing-down candidate blackleg resistance genes. This study provides a novel genomic resource towards the identification of candidate genes for breeding disease resistance in B. rapa and its relatives.